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ABSTRACT 

We use a dynamical systems approach to study thawing quintessence models, using a multi- 
parameter extension of the exponential potential which can approximate the form of typical 
thawing potentials. We impose observational constraints using a compilation of current data, 
and forecast the tightening of constraints expected from future dark energy surveys, as well 
as discussing the relation of our results to analytical constraints already in the literature. 
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1 INTRODUCTION 

Thawing quintessence models, in the terminology of Caldwell & 
Linder (2005), are those in which at early times the dark energy 
has a much lower density than matter, but where the dark energy 
density begins to evolve once it becomes a significant fraction of 
the total (Steinhardt, Wang & Zlatev 1999; Huey & Lidsey 2001; 
Nunes & Copeland 2002). It is characterized by an equation of state 
which is initially w = — 1 to high accuracy, and which then 'thaws' 
to w > —1 by the present. Provided the thawing process is suffi- 
ciently slow, these models are in agreement with current observa- 
tional data. 

At first sight the initial conditions for thawing quintessence 
appear very unnatural, as they require the same fine-tuning of a 
small initial density that one finds with the cosmological constant. 
This fine-tuning could presumably be explained using the same 
sort of anthropic string landscape argument often invoked for a 
pure cosmological constant (Bousso & Polchinski 2000; Susskind 
2003). However the case of thawing quintessence may in fact be 
more appealing, because at least for very steep potentials one can 
argue that quantum fluctuations of the quintessence field acquired 
during early Universe inflation may drive the field to low energy 
densities (Malquarti & Liddle 2002). A complete cosmology from 
inflation through to the present may well lead to a significant prob- 
ability of thawing quintessence behaviour, without imposing addi- 
tional anthropic constraints. 

In light of this, it is useful to characterize the types of thaw- 
ing model allowed by present data, and indeed there have been 
several papers recently exploring aspects of thawing quintessence 
models, though these have largely ignored the issue of natural- 
ness of initial conditions. Inflationary models can be character- 
ized by the slow-roll approximation, but as this is not generally 
valid for quintessence (Bludman 2004; Capone, Rubano & Scud- 
ellaro 2006; Linder 2006; Cahn, de Putter & Linder 2008) some 
papers have sought an analogue to this approach for dark energy 
l |Scherrer& Sen 20081 IDutta & Scherrer 2008t . Another method, 
employed by Crittenden, Majerotto & Piazza (2007), involves 
smoothness requirements on the potential, while other work has 



developed a classification of dark energy models by means of 
a calibration of their time variation (Linder 2006; Linder 2008; 
ICahn et al. 20081 Ide Putter & Linder 20081 . Analytical bounds on 
thawing potentials have previously been derived by Caldwell & 
Linder (2005), Linder (2006) and Barger, Guarnaccia & Marfa- 
tia (2006), who also considered previous observational constraints. 
Expansions of the putative quintessence potential, either using the 
flow equations (Huterer & Peiris 2007) or directly (Sahlen, Lid- 
dle & Parkinson 2007) have been used to explore broad classes of 
models to contrast the thawing and non-thawing regimes. 

In this paper we analyze thawing models from a somewhat 
different perspective. We look at a multi-parameter potential fam- 
ily V(tf>), which includes the exponential potential as a special case, 
and which has enough freedom that the family can represent arbi- 
trary values of V, dV/d(f>, and d 2 V/d(f> 2 at the initial field value. 
Reminiscent of the slow-roll approximation to inflationary observ- 
ables, which depend only on the potential and its first few deriva- 
tives, we aim to capture with this potential the full spectrum of 
thawing quintessence phenomenology. This is then confronted with 
present observational data, and the capabilities of future data as- 
sessed, before a consideration of the analytical constraints derived 
in previous analyses of thawing behaviour. 



2 THAWING QUINTESSENCE 
2.1 Quintessence dynamics 

We assume a minimally-coupled scalar field (j> acts as quintessence. 
The pressure/density relation is w = p/p with 



(i) 



where dots denote derivatives with respect to time and V(<fi) is 
the field's self-interaction potential. The variation in <j> obeys the 
Klein-Gordon equation 



< + 3H<j> + V' = 0, 



(2) 
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where the prime indicates a derivative with respect to <f>. This shows 
that the field rolls down the slope of its potential V{4>) and that its 
motion is damped by the Hubble parameter which, in units with 
8nG = 1, may be written 



(3) 



where a is the scale factor of the Universe and pt is the total density. 
We assume a flat matter-dominated Universe throughout. 

The following dynamical systems approach, equations {4} 
through fllOt , was developed by Copeland, Liddle & Wands (1998), 
de la Macorra & Piccinelli (2000), Ng, Nunes & Rosati (2001), and 
Scherrer & Sen (2008). Equations ([2} and I0 may be rewritten in 
terms of three new variables; x, y and A, defined as 



y/6H' 



V(4>) 

3H 2 ' 



(4) 



For a Universe containing only matter and a scalar field, the system 
may be written as ( Ng et al. 200T) 



aUv = -3^ + A^ 2 + ^(l + x 2 -y 2 ), 
dy 3 3 . 2 2n 



(5) 

(6) 
(7) 



where N is the logarithm of the scale factor, N = In a, and V is 
defined as dSteinhardt et al. 1999) 



VV" 

iY'Y 



(8) 



For typical potentials, including the one we introduce below, V is 
a function of A (though only expressible as such if \(<f>) has an 
analytic inversion) and hence the system is an autonomous one (see 
e.g. Fang et al. 2008), though this property is not needed in our 
analysis. 

In a Universe with flat geometry, the density parameter of the 
scalar field becomes 



2 . 2 

x +y , 



and its effective equation of state is 



7 = 1 + w — 



2x 2 



x 2 + y 2 



(9) 



(10) 



2.2 Generalized thawing potential 

The initial value of the scalar field is assumed to be fixed at some 
early time in the matter-dominated era, and can be set to zero with- 
out any loss of generality, ie. cj>i = 0. Having made this choice, 
we can then explore a wide range of thawing models by taking the 
quintessence potential to be of the form 



V{<t>) = Vi exp(-c0)[l + a<? 



(11) 



where Vi is the initial value of the potential, and c and a are pa- 
rameters. This form is useful because of two properties. Firstly, it 
reduces to the well-known exponential potential in the case a — 0; 
this property motivates our choice in favour of a simple Taylor ex- 
pansion around the origin, as often used in direct reconstruction of 
the quintessence potential (e.g. Sahlen, Liddle & Parkinson 2005). 
Secondly, in the vicinity of the origin it can approximate arbitrary 
potentials up to their second derivative. In the limit a — *• c it can 




Figure 1. The evolution for Ai = 0.1 (outer curve) to 1 (inner curve) in 
steps of 0.1 for c = 4. The dashed semicircle represents the present-day 
dark energy density parameter Q.^ « 0.75. Trajectories which reach this 
may correspond to the evolution of the Universe. The origin is the initial 
point at some early time during matter domination, and the unit semicircle 
represents domination by the scalar field. 



also approximate the hilltop quintessence models recently studied 
by Dutta & Scherrer (2008). 

Our results will therefore have two separate interpretations. 
The first is to consider equation Jilt to be an exact potential (some- 
what similar to the Albrecht-Skordis (2000) potential), valid for all 
<f), enabling generalization of results from the exponential case. The 
second is to consider the results to be valid for all thawing poten- 
tials, provided the variation in <f> is small enough that our potential 
is a good local approximation to an arbitrary potential. We will take 
the convergence condition for this potential to be \atf>\ < 0.1. We 
also restrict the parameters to c > a so that initially V' < and 
hence the field value increases with time. 

Putting this potential into equation (|7) leads to 



and the initial value of A is given by 



Ai 



(12) 



(13) 



We now allow the system to evolve out of an initial perturbation 
from cj>i = 0, exploring the parameter space by plotting the evolu- 
tion of the system over time for a range of values of a and c (e.g. 
see Fig. [I). Equation $9^ implies that < x 2 + y 2 < 1 and so the 
resulting trajectories are confined to within the unit circle. From 
the definitions of x and y, equation {4}, it can also be seen that 
making y negative and reversing the temporal coordinate (and thus 
— * — cj}), results in symmetrical behaviour below the axis to that 
above. Only the upper half-disc therefore need be considered for 
a full picture of possible trajectories. Furthermore, the lower half- 
disc corresponds to a contracting Universe, which is ruled out by 
observations as far as the past and present Universe is concerned. 

Our initial condition keeps the initial scalar field velocity at 
zero, which is not quite appropriate. Cahn et al. (2008) carried out 
an analysis of the early-time behaviour of quintessence fields to 
compute the leading-order early-time behaviour of the velocity. In 
our case this translates as x\ oc yf, meaning that the initial veloc- 
ity is indeed highly suppressed. We confirmed numerically that the 
difference to the trajectories is completely negligible. 
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Figure 2. Families of potentials which reach = 0.75, and have positive 
values of the parameter a, for c = 0.5 (top solid line) to 5 (lowest solid 
line) in steps of 0.5. The short-close-dashed line delineates the range of val- 
ues allowed by recent observational data and the dotted line represents the 
limit where a = 0. The dot-dash line indicates the ten-year observational 
prospects discussed in Section 13.21 while the spaced-dashed and the long- 
dashed curves represent potentials with T; = and To = respectively. 



3 OBSERVATIONAL CONSTRAINTS 
3.1 Current constraints 

The present-day dark energy density parameter has an observed 
value of ~ 0.75 (Komatsu et al. 2008) and using equation ([9} 
we can determine which trajectories reach this value and as such 
may be representative of the real world (e.g. Fig.[T|l. 

The set of possible models is better constrained, how- 
ever, by considering their representation in the Wo-w a plane. 
Thawing dark energy is well parameterized by the equation 
dChevallier & Polarski 200TllLinder 20031 



w(a) — wo + w a (l — a), 



(14) 



where wo represents the present-day value of w, and w a determines 
the change in w with the scale factor a. Equation dlOt can now be 
used to derive expressions for wq and w a in terms of the variables 
of the autonomous system. Since a = 1 at the present day we have 



w = 



(15) 



where xo and yo are the present-day values of x and y. The obser- 
vational constraint that w is close to —1 therefore requires viable 
trajectories to remain close to the j/-axis at the time they reach the 
semicircle indicating the correct dark energy density. 

The time dependence of the equation of state is also obser- 
vationally constrained, though less strongly. Combining equations 
dlOt and d 1 4b leads to the expression for w a 



Wa 



4:X y 
\xl + ylf 



-?>x y + A 



3 2 



(16) 

where An is the present-day value of A. 

The coordinates at which trajectories reach the present dark 
energy density can now be used in these two equations to plot the 



Figure 3. As Fig. [2] but for potentials with negative values of the parameter 
a and c = (lowest solid curve) to 1 (top solid line) in steps of 0.1. The 
spaced-dashed line is the — w a > 1 + Wq thawing limit, while the long- 
dashed curve represents the curve where X = 3/4 (see Sections- 



predictions from families of potentials in the wo-w a plane, shown 
in Figs. |2] and [5] These values are compared to the observationally- 
allowed region (short-close-dashed curve) obtained by Komatsu et 
al. (2008) from a combination of baryon acoustic oscillation, type 
1A supernovae and WMAP5 data at a 95% confidence limit, and 
would only be marginally improved by the inclusion of big-bang 
nucleosynthesis constraints. 

One might further worry whether the approximation of con- 
stant w a is accurate enough, and indeed Dutta & Scherrer (2008) 
indicated that in some parameter regions it will not be. We quantify 
this below by delineating the regions of parameter space in which 
w a has varied by more than 25% from redshift one to the present. 

Each of the solid curves in Figs. [2] and [3] represents a family 
of potentials with the same value of the parameter c in the expo- 
nent, but different values of the linear parameter a and therefore 
Ai (remember Ai = c — a). They are not the evolutionary tracks 
of particular potentials, but rather indicate the points in the wo-w a 
plane which different potentials have reached at the present day. 
The largest values of a are at the left-hand end of the curves in 
Fig. [2] where they all share a common origin at the point [—1,0], 
corresponding to the cosmological constant. Potentials at this point 
have c — a which means that Ai = and so they are flat, while 
nearby potentials have negative curvature, ie. F < 0, and corre- 
spond to the case of hilltop quintessence as described by Dutta & 
Scherrer (2008). Moving to the right along the curves the potentials 
have decreasing values of a and hence increasing Ai, ie. increas- 
ingly negative initial slope. 

Useful orientation into the behaviour of the family of models 
comes from studying the time dependence of T. We see from equa- 
tion J 12b that A is monotonically increasing with time, and hence F 
too is an increasing function whose initial value may be positive or 
negative (we keep c > throughout, corresponding to cj> increasing 
with time). The sign of F is the same as of V" , and hence deter- 
mines whether the potential is steepening or become more shallow 
with time, with V — corresponding to a point of inflection in the 
potential. We can then classify the models as follows. If the ini- 
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tial value of F, Fj, is already positive, then Y stays positive forever 
and the potential is becoming shallower with time (the exponen- 
tial case is the archetypal example). If the present value of Y, To, 
is negative, then Y has been negative throughout the past evolu- 
tion of the Universe, corresponding to a steepening potential as in 
hilltop quintessence. If neither of these conditions is satisfied, then 
r changes from negative to positive during the evolution, with the 
potential initially steepening and later becoming more shallow. 

In the figures, the curvature of the potentials increases to the 
right and eventually potentials with To = (long-dashed) and 
subsequently 1^ = (spaced-dashed line) are reached, corre- 
sponding to points of inflection in the potential curve. After the 
points of inflection the slope of the potentials increases along the 
curves, which then all head towards the exponential limit as a ap- 
proaches zero. They eventually graze the exponential curve where 
c — \ — constant and the curvature reaches its maximum value of 
Yi = 1. Higher values of c than those plotted in Fig.[2]are possible, 
but only as long as c ~ a. For example, the trajectory of c = 20 
in the x-y plane only reaches the dark energy density parameter 
curve when c — a < 10~ 10 . 

The curves in Fig.[3]represent similar families of potentials to 
those in Fig.|2]but with negative values of the parameter a. These 
curves begin from the exponential limit at their left-hand end and 
have increasingly negative values of a to the right. This means 
that they have increasing Ai, ie. the initial slope of the potentials 
becomes ever steeper to the right. Here the curves once again di- 
verge from the exponential case as the magnitude of a increases. 
This is a consequence of their increasing linearity, which causes 
the curves to asymptotically approach the family of linear poten- 
tials with c = r = 0, as eventually Ai ~ c. Eventually each family 
of curves comes into conflict with observations, which therefore 
give a lower limit on a for each value of c. The region of poten- 
tials in the c-a parameter space consistent with the observational 
constraints is shown in Fig. [4] 

The range of potentials for which convergence remains good 
up until the present day can be found by rewriting the convergence 
condition as 

0.9 < — ^— < 1.1, (17) 

c - A 

and is also shown in Fig. [4] Potentials with c < are not consid- 
ered here, since they simply correspond to those for which the field 
rolls in the opposite direction and so add nothing new to the anal- 
ysis. Figure 4 also shows the region (essentially a > 1) where the 
assumption of constant w a used in our observational comparison is 
breaking down; the observational results should be considered less 
robust in this region. 

Admitting potentials with a < allows V ((f)) to become neg- 
ative in the future, which may prompt future recollapse of the Uni- 
verse (Kallosh et al. 2003). Such potentials are included in our anal- 
ysis, though they may be considered as being approximate repre- 
sentations of the behaviour of arbitrary potentials from early times 
to the present day, without necessarily extrapolating into the nega- 
tive energy density regime. 

The representation of our results in terms of the more familiar 
slow-roll parameters from inflation, e and 77, is easily found from 
their definitions (in natural units) 

- - m'- 1 **- 

Th = ^r=2cAi-c 2 , (19) 




a 



Figure 4. The allowed potentials in the c-a parameter space. The region 
beneath the short-dashed line is allowed by the observational constraints, 
while the spaced-dash line below it gives an approximation to the ten- 
year observational prospects discussed in Section fOl The region between 
the dotted lines contains potentials for which convergence remains good 
(\arp\ < 0.1) up until the present day and the solid line represents the phys- 
ical limit for which A is positive (i.e. dV/d(f> remains negative). Models to 
the left of the dot-dashed line have w a constant to a good approximation; 
our observational treatment is less reliable for those to the right of this line. 
Finally, the long-dashed line is where r\ = 0, which divides the parame- 
ter space into shallowing potentials on the left and those which undergo a 
period of steepening on the right. 



and is shown in Fig. [5] As would be expected from the results of 
Scherrer & Sen (2008) and Dutta & Scherrer (2008), the assump- 
tion of constant w a fails for large negative 77 (the top panel), mean- 
ing our observational comparison method is not robust, but is fine 
otherwise. 



3.2 Future prospects 

As well as the present-day observational constraints, we can use 
predictions of future bounds to give a sense of how the range of pos- 
sible potentials may be further restricted. Each graph (except Fig.[6jl 
displays an approximation to the ten-year observational prospects 
(dot-dashed curve) under the assumption that ACDM is the true 
model. This was achieved by approximating the predicted future 
Wo-Wa limits from the literature and then resizing the present ob- 
servational constraints curve to fit them, thereby assuming similarly 
distributed future constraints. These constraints roughly correspond 
to predictions for both forthcoming Stage 3 programmes, such as 
the Dark Energy Survey, at 68% c.l. <Abbott et al. 2005) and future 
Stage 4 projects, such as JDEM, at 95% c.l. JAlbrechtet al. 2006). 
Figs.|2]and|3]show that a significant limiting of the allowed region 
in the wo-w a plane is possible over the next ten years, but Fig. [4] 
shows that this does not translate into a correspondingly large re- 
duction in the c-a parameter space used in our model. There is 
however a slightly greater reduction in the size of the allowed re- 
gion in the rj-e plane, as shown in Fig. [5] 
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Figure 5. The allowed region of potentials in the rf—e plane. Positive a with 
negative I\ is shown on at the top, positive a with positive T,; is shown in 
the middle, and negative a (all of which have V{ > 0) are at the bottom. 
The lower two figures correspond to potentials which have been shallowing 
throughout their history. The region below the dashed line is consistent with 
present observations and the space-dash line gives an approximation to the 
ten-year observational prospects discussed in Section |X2l Potentials below 
the dotted line all have \ot<j>\ < 0.1 up until the present day, and so converge 
well. The e axis is where T; = and the rj axis is where a = c, while the 
solid line in the lower two figures represents the limit where a = 0. In the 
top figure, only models to the right of the dot-dash line have w a constant 
to a good approximation (this condition is always satisfied in the lower two 
graphs). 




0.80 



Figure 6. The solid line is the upper thawing bound of — w a > 1 + Wo, 
the long-dashed line is the Barger et al. (2006) limit and the short-dashed 
curve is where X = 3/4. The dot-dashed line is the original Caldwell & 
Linder (2005) lower thawing bound and the dotted line is the curve where 
X = 3/2. 



4 ANALYTICAL CONSTRAINTS 

A number of different analytical constraints have been placed on 
thawing models and the present situation could perhaps bear some 
clarification. Fig.|6]shows constraints in the wo-w a plane from pre- 
vious analyses of thawing dynamics. The upper limit holds in all 
cases, as only below this line is the requirement w > — 1 fulfilled, 
ie. — w a > 1 + Wo (solid line). Also, from the original definition of 
thawing behaviour by Caldwell & Linder (2005), comes the lowest 
of the linear boundaries, — w a < 3(1 + Wq) (dot-dashed line). This 
comes from the requirement that <f>t < <j} f° r potentials where <f> is 
decreasing for all time. The slope of such potentials is monoton- 
ically increasing (becoming less negative) as the field asymptoti- 
cally approaches the minimum of its potential curve. 

Linder (2006) defined two useful ratios of the Klein-Gordon 
terms and the relationship between them 



X 



H4> 



Y 



V 



X 



Y 



l + Y 



(20) 



The Klein-Gordon equation can be used along with these defini- 
tions to obtain an expression for X in terms of x, y and A 



X 



^6; 



(21) 



This allows us to plot given values of X using the generalized po- 
tential. The aforementioned lower bound of — w a < 3(1 + wo) 
was shown by Linder (2006) to be an approximation to the more 
precise limit of A = 3/2 (dotted curve), although this curve does 
still contain the implicit approximation that H w 2/3t. Barger et 
al. (2006) applied the original Caldwell & Linder (2005) bound at a 
redshift of 1 , allowing for a significant tightening of the constraints 
on potentials with monotonically increasing slope and showed that 
such potentials satisfy — w a > (3/2)(l + Wo) (long-dashed line). 
This is again a linear approximation and we plot this bound as the 
more accurate curve where X — 3/4 (short-dashed curve). 

Potentials which steepen as the field rolls down may still sat- 
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isfy the basic requirement for thawing that w begins very close 
to —1 before increasing at later times, despite permitting a more 
complex thawing behaviour. In principle their range extends all the 
way to the wq — — 1 line, which represents the divide from the 
present-day phantom regime but, as mentioned in Section [JT1 in 
our model such potentials need to be increasingly finely tuned the 
closer they are to wo = — 1. Linder (2006) mentioned that po- 
tentials which steepen as the field rolls down should lie below the 
bound of <j)t < cj>. We find that this requirement however, whilst 
sufficient to ensure that potentials steepen as the field rolls down, 
is not a necessary condition. We find that potentials which do un- 
dergo a period of steepening slope since the onset of matter dom- 
ination all lie beneath the curve where T\ — (spaced-dashed 
line) in Fig. [2] This represents the limit below which the potentials 
have all reached the point of inflection in the generalized poten- 
tial curve by the initial early time during matter domination. We 
find that this coincidentally corresponds to a very good accuracy to 
the aforementioned Barger et al. (2006) linear bound, which may 
therefore be used as an upper limit for potentials which have not 
had monotonically-increasing slope throughout matter domination. 

The long-dashed curve in Fig. [2] shows potentials for which 
To — 0, ie. they reach their point of inflection at the present day, 
and as such have been steepening throughout matter domination. 
Potentials between this curve and the Ti = line are therefore 
those which have crossed their point of inflection at some time dur- 
ing the matter-dominated epoch, and so the distance between these 
two curves is related to the age of the Universe. It should also be 
noted that above the X — 3/4 curve (long-dashed in Fig. [3j there 
are always two potentials at each point in the wo-w a plane, and 
above the Ti = curve (spaced-dashed in Fig.O there are always 
two potentials with monotonically increasing slope at each point. 

Potentials above the Ti = curve in Fig. [2] and those with 
negative a (ie. all those in Fig. [3}, have all had positive curvature 
(Ti > 0), and as such increasingly shallower slope, since the ini- 
tial time. These are the potentials to which the Caldwell & Linder 
(2005), Linder (2006), and Barger et al. (2006) limits all apply. 
We find that the lower limit on such potentials is the case where 
c = 0, i.e. the family of linear potentials. Since linear potentials 
lie approximately along the curve of X = 3/4 we find that this 
could therefore be used as a lower bound on potentials with mono- 
tonically increasing slope to which the previous limits have been 
applied, as shown in Fig. [3] (long-dashed curve). Also shown in 
Fig. [3] is our finding that the exponential potentials (large-dotted 
curve) represent an upper limit and may still be best approximated 
by the original Caldwell & Linder (2005) upper bound (spaced- 
dashed line). 



5 CONCLUSIONS 

The form of the scalar field self-interaction potential studied 
here allows observational constraints to be placed on thawing 
quintessence models, encompassing a wide range of thawing be- 
haviour. A pure exponential potential for example, which has the 
greatest observationally-allowed slope (Ai ~ 1), is currently re- 
stricted to having c less than about 1, and upcoming observations 
could reduce this by up to a third. Despite applying strict con- 
vergence conditions, we expect the chosen potential to be fairly 
reliable within the whole of the observationally-allowed region. 
This is especially true when considering the ten-year observational 
prospects, for which the value of \ot<f)\ peaks at around 0.3. 

The analytical constraints on thawing models in the wo-w a 



plane are tightest for potentials which have had monotonically in- 
creasing slope throughout their evolution history. Their range may 
be reasonably well approximated using limits from previous analy- 
ses, but within these bounds there are always two potentials at any 
given point in the wo-w a plane. Potentials which have undergone a 
period of decreasing slope have an upper bound which is lower than 
that for purely shallowing potentials, and which is coincidentally 
well approximated by an existing limit, but only fine-tuning argu- 
ments restrict these potentials otherwise. Future data will require 
new analytical constraints to be fitted to potentials in the thawing 
region for the sake of accuracy, but for now existing limits properly 
applied offer a fairly accurate representation of the true range of 
viable models. 
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